Infrared divergences, mass shell singularities and gauge dependence of the dynamical 

fermion mass 
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We study the behavior of the dynamical fermion mass when infrared divergences and mass shell 
singularities are present in a gauge theory. In particular, in the massive Schwinger model in covariant 
gauges we find that the pole of the fermion propagator is divergent and gauge dependent at one 
loop, but the leading singularities cancel in the quenched rainbow approximation. On the other 
hand, in physical gauges, we find that the dynamical fermion mass is finite and gauge independent 
at least up to one loop. 
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The pole of the propagator for a particle defines the 
mass of the particle and, in a gauge theory, the pole is 
expected to be gauge independent P, Intuitively it 
is clear that the physical mass of a particle should be 
gauge independent and quantitatively it can be argued 
as follows. Let us consider a fermion of tree level mass m 
interacting with a photon field in a covariant gauge. Then 
the physical mass m p of the fermion is determined from 
the vanishing of the complete two point function S^ 1 = $ 
— m — S(p, m) = 0. Here E(p, m) denotes the fermion 
self-energy to all orders. At one loop, for example, the 
fermion self-energy (see Fig. [TJ in /^-dimensions leads to 
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FIG. 1: Fermion self-energy in QED at one loop, 
the equation 
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where £ denotes the (covariant) gauge fixing parameter 
and 
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Since to one loop order = r 

right hand side of (fTJ) is already of order e 2 , we can let 
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in the terms — m) on the right hand side 
which shows that the right hand side vanishes at the pole 
j> = vfip if the functions F\,F% are well behaved. Fur- 
thermore, if we can Taylor expand the self-energy around 
the pole to write 

£« (p, m) = (m W - m) + - m«) (£« )' + • • • , 

p—rrtp ' 

(3) 

only the first term survives at the pole. As a result, at the 

9m (1) 

pole |T]) leads to g| = showing that the pole of the 
propagator (mass of the particle) is gauge independent. 
Gauge independence of the pole at higher loops is best 
studied through an analysis of the Nielsen identity for 
the two point function which we will discuss later [3j-l5( . 

It is clear from ([2]) that the functions F\ , F2 become 
infrared singular for space-time dimensions D < 3 so 
that the right hand side of (fT]) may not vanish in gen- 
eral. In addition, as we will see in lower dimensional 
models, mass shell singularities develop in the fermion 
self-energy so that a Taylor expansion (as in ©) is not 
possible and the gauge independence of the pole of the 
propagator can not be guaranteed. These are two possi- 
ble difficulties that can invalidate the formal argument of 
gauge independence given above and when infrared diver- 
gences and/or mass singularities are present, the question 
of gauge independence of the pole has to be examined 
anew. This is important since the effective action (and 
its imaginary part) for fermions interacting with a gauge 
field background depends on the fermion mass [6j and 
being a physical quantity (the imaginary part is related 
to the vacuum decay rate), it cannot possibly be gauge 
dependent. 

Infrared singularities are well known to be a problem in 
lower dimensions 0, but even in 4-dimensional QCD, 
they can arise (in the fermion self-energy) at two loops 
and beyond because of gluon loops @ . Of course, intu- 
itively it is clear that a strong infrared divergence signifies 
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possible confinement of the fermion in the theory with- 
out any physical fermion asymptotic states. (As we know 
this is true even in Abelian theories in lower dimensions 
where the strong infrared behavior of the Coulomb force 
leads to confinement of fcrmions SEMI-) Nonethe- 
less one can calculate the fermion mass for the perturba- 
tive asymptotic states as is done in the case of "current" 
quark masses which are then used in various calculations 
and it is worth investigating whether such a calculation 
can lead to a gauge independent result. 

For these reasons we have chosen to study this ques- 
tion systematically. We find that in covariant gauges 
infrared divergences as well as mass shell singularities de- 
velop in the fermion self-energy at one loop when D < 3 
invalidating the standard arguments for gauge indepen- 
dence of the pole of the propagator. However, we show 
that the leading divergences cancel when the fermion self- 
energy is summed to all orders in the quenched rainbow 
approximation [l4], EH . This signals the fact that the in- 
frared divergences and mass shell singularities may not be 
present in a physical gauge such as the axial and Coulomb 
gauges. We carry out such a calculation for the 1 + 1 di- 
mensional massive QED and show explicitly that there 
are no infrared divergences in the self-energy up to two 
loops. We also find explicitly that the pole of the propa- 
gator is gauge independent at least up to one loop. 

The one loop fermion self-energy in Fig. [T] can be cal- 
culated in a general covariant gauge in a straightforward 
manner and in dimensional regularization it has a closed 
form expression in terms of hypergeometric functions. 
However, to see the singularity structure, we note that 
near the mass shell p 2 — > m? it behaves like 



This is infrared divergent, has a mass shell singularity 
and is also explicitly gauge dependent. Although the 
leading divergence in ([5]) can be set to zero with the 
choice £ = — 1, subleading divergences (which correspond 
to mass shell singularities) persist. In fact, there is no 
value of £ which will make the self-energy completely 
well behaved. It is rather counterintuitive that the mas- 
sive Schwinger model has an infrared divergence while 
the massless Schwinger model does not. This has its ori- 
gin in the fact that the Dirac algebra in the self-energy 
in Fig. [T] leads to the numerator (say, in the Feynman 
gauge for illustration purposes) 7 M (^ + m)j^ = (— 2e^ 
+ (2 + 2e)m). When m ^ 0, this numerator is finite in 
the limit e — > and does not moderate the infrared di- 
vergence arising from the integral while when m = 0, the 
numerator behaves as —1tj> and the e in the numerator 
moderates the infrared behavior making it finite. 

Since the one loop self-energy is divergent, it is worth 
asking if the higher loop corrections can cure the prob- 
lem. A general analysis of this question is rather difficult. 
Instead we will show that when summed to all loops in 
the quenched rainbow approximation, the leading diver- 
gence in the fermion self-energy does indeed vanish. First 
let us recall that the overlapping divergence graphs, in 
general, have a softer infrared behavior [l6| so that in 
studying the leading divergence structure these can be 
neglected. Second, let us note here that when infrared 
divergences are present, normally including the polariza- 
tion tensor in the photon propagator at higher loops im- 
proves the behavior. However, in the case of the massive 
Schwinger model, the one loop polarization tensor has 
the explicit form 
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This shows that for D > 3, there are the usual ultraviolet 
divergences (coming from T(3 — D) and T(2 — §)), but 
there is no mass shell singularity so that the self-energy 
can have a proper Taylor expansion as in ([3]). On the 
other hand, when D < 3 there are not only infrared diver- 
gences (the gamma function singularities in these cases 
can be identified with infrared divergences) but also mass 
shell singularities which show that a Taylor expansion of 
the self-energy is not possible. Explicitly, for example, we 
see from ^ that for D — 2 + 2e the leading divergence 
of the self-energy has the form 
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It is clear from that when the fermion mass vanishes 
Hr(p 2 ) = — — giving a mass to the photon in the (mass- 
less) Schwinger model. When m ^ 0, however, the func- 
tion Hr(p 2 ) is non-analytic. In fact, for < p 2 < Am 2 
we obtain from ([7]) 
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while for p 2 > Am 2 we have from {7J 
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Equation © allows us to take the massless limit yielding 
the result for the Schwinger model. On the other hand, 
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for a nonzero fermion mass, it is 
take the limit p 2 — > leading to 



which allows us to 



n T (p 2 )^>--[i-i] = o, (to) 

showing that a photon mass is not generated. (The ab- 
sence of a "Higgs" phenomena in the massive Schwinger 
model is already known from a different analysis [l^ . [l3| .) 
As a result, including vacuum polarization in the photon 
propagator in higher orders does not help and we restrict 
our analysis to the quenched rainbow approximation (see 
Fig. EJ. 




FIG. 2: Fermion self-energy diagrams in the quenched rain- 
bow approximation. 

The leading divergence terms in the self-energy dia- 
grams in Fig. [2] for the massive Schwinger model can be 
calculated in a straightforward manner and the contri- 
bution at the (n + l)th loop turns out to be recursively 
related to the contribution at the nth loop in a simple 
manner 
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This simple recursive behavior of the leading divergences 
in the quenched rainbow approximation allows us to sum 
the self-energy to all orders (in this approximation) and 
yields 
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This shows explicitly that the leading divergence van- 
ishes in the self-energy when summed to all orders in the 
quenched rainbow approximation (for any value of £). 
This simple relation can also be seen in another way. We 
note that the leading term in the self-energy to all orders 
can be written as an integral equation ("gap" equation) 
of the form 



X L (p) = ^(p) 
d D k .. 
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The leading infrared divergent contribution comes from 
the region k — > and as j> — > m this equation is solved 
in this regime by = — — m) —> 0. (We have 



explicitly seen the same behavior in scalar QED in 1 + 1 
dimensions in covariant gauges as well.) 

Of course, this does not say anything about the sub- 
leading divergences which, in principle, are much harder 
to calculate since they involve the overlapping graphs. 
However, the vanishing of the leading divergence (when 
summed) suggests that the fermion self-energy may be 
better behaved when calculated in a physical gauge and 
this is what we discuss next. This can be understood as 
follows. In the covariant gauge, the fermion self-energy 
receives contribution from the unphysical as well as the 
physical degrees of the photon and a resummation may 
be necessary to see the cancellation between the contri- 
butions from these degrees of freedom. However, having 
seen the cancellation in the covariant gauge, it is clear 
that since a physical gauge has only physical degrees of 
freedom, there is no room for cancellation any more and 
the result should already be well behaved at one loop. In 
fact physical gauges, such as the axial gauge, are known 
to tame the infrared divergence much better and if there 
is no infrared divergence divergence at one loop (as we 
will see shortly), summing the rainbow graphs would no 
longer be a meaningful approximation. 

Let us consider the massive Schwinger model in the 
homogeneous axial gauge n ■ A = where n M denotes a 
fixed direction in space-time. In this case, the photon 
propagator has the form 
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and satisfies 
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The poles of the propagator in (|T4j) at n ■ p = are 
conventionally handled through the PV (principal value) 
prescription [171 1 1 81 ] which is known to soften the infrared 
behavior. It is also known that the PV prescription may 
be problematic for time-like (n 2 > 0) and light-like 
(n 2 = 0) axial gauges. For this reason we restrict our- 
selves to the class of space-like (n 2 < 0) axial gauges. It 
is worth noting here that in D = 2, the space-like ax- 
ial gauges also include the Coulomb gauge. This can be 
seen in two different (but equivalent) ways. First, we 
note that if we choose the gauge A\ = or n M = (0, 1), 
the only nontrivial component of the propagator in (|14[) 
is given by 



«A)o(p) 



2 ' 

Pi 



(16) 



which coincides with the propagator in the Coulomb 
gauge. Alternatively, we note that even after we im- 
pose the Coulomb gauge d\A\ — which implies that 
A\ = Ai(t), there is a residual gauge symmetry preserv- 
ing this gauge choice, namely, 



A 1 ^A[=A 1 (t) + d 1 A(t,x), 9i^i=0. 



(17) 
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This determines A(t,x) = xa(t) + f3(t) where a(t),(3(t) 
are arbitrary functions of t. This leads to (see (fT7|0 

A' 1 = A 1 {t) + a{t), (18) 

and shows that if we choose a(t) = —Ai(t), A\ = which 
corresponds to the axial gauge. Coulomb gauge is known 
to be notoriously singular in D = 2. However, in con- 
sidering this as belonging to the class of axial gauges in 
D = 2, it inherits the PV prescription which gives it an 
improved behavior. 



In D = 2 the two vectors 

n>*,n%, = ei M 'n v , (19) 

define a complete basis and, therefore, any vector (includ- 
ing the Dirac matrices 7^) can be expanded in this basis. 
Using this the calculation of the fermion self-energy at 
one loop (see Fig. [1} can be carried out in a manifestly 
covariant manner. The PV prescription indeed makes 
the integral infrared finite and we obtain 



Y^^(p, m, n) = h 
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where we have defined 



\J in ■ p) 2 — m 2 n 2 

The seemingly divergent second term in (|20[) . in fact, van- 
ishes as n ■ p — > so that the self-energy is completely 
well behaved. We also note that the complete gauge de- 
pendence (n M dependence) of the self-energy is contained 
in the second term which vanishes as j) — ► m so that 
the one loop correction to the pole of the fermion prop- 
agator is gauge independent. This calculation can also 
be extended to two loops. The axial gauge calculations 
are, in general, extremely difficult at two loops and be- 
yond. We have not yet been able to evaluate the two loop 
self-energy explicitly. However, analyzing the integrand 
at this order (including overlapping graphs) we find that 
the integral appears to be infrared finite. But, without 
an explicit evaluation of the self-energy at two loops, it 
is not yet clear whether the pole of the propagator will 
be gauge independent at this order as well. 

As we have mentioned earlier, the gauge dependence of 
the pole of the propagator is best analyzed using Nielsen 
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FIG. 3: The two three point vertices F M , G M . 



identities which can be derived using the BRST symme- 
try [2(3] of the theory. Let us briefly sketch here how it is 
derived in a homogeneous axial gauge. The gauge fixed 
QED Lagrangian density (including free ghosts) in any 
dimension is given by 

C = -\f^F^ + ij ( 7 "(i5 /1 - eA^) - m) V 

+ Fn ■ A + cn ■ dc, (22) 

where F denotes an auxiliary field introduced for gauge 
fixing. We can introduce the Lagrangian density of the 
standard sources needed to study BRST (Ward) identi- 
ties. However, in order to study the gauge dependence 
of the fermion two point function we need to introduce 
an additional source (fermionic with ghost quantum 
number) for the composite field (cA^). Thus the com- 
plete Lagrangian density for the sources is given by 

Source = J^Ap + JF + i(i>x ~ XVO + *(c?y - fjc) 

+ ie{Mcip-$cM) + H f *cA fll (23) 



/ 
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where M,M denote respectively the sources (bosonic spinors with anti-ghost/ghost quantum numbers) for the 
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composite BRST variations of the fields V> and ip. With 
these it is straightforward to derive the identity describ- 
ing the gauge dependence of the fermion two point func- 
tion 



J _af3 



F^a^p, -p)S^J{p)+S a ^(p)G f _ inl3 {Q, -p,p), 

(24) 

where a, /3, 7 denote the Dirac indices and F^, G M denote 
the three point vertices shown in Fig. [3] Using the form 
of the photon propagator in (Ti"4"|) . it is straightforward to 
check that 



n ■ p 



D X „(P) 



Pu 

n ■ p 



DM- 



(25) 



Furthermore, using this relation (Eq. ([25j)) we can ex- 
plicitly check graphically up to two loops that the Nielsen 
identity (|2"41 holds which is consistent with our earlier ob- 
servation that there is no infrared divergence at least up 
to two loops. 

The Nielsen identity (|24|) can be studied order by or- 
der in a loop expansion. We note from (j2"3")l that there 
are no tree level three point vertices corresponding to 
Ffj, >a r/,Gfj, t ar/ (see Fig. [3]). They can arise only through 
radiative corrections and this is consistent with our ex- 
pectation that gauge dependence in the fermion mass can 
possibly arise at the loop level. The vertices F M , Q7 , G>, Q7 
indeed arise beginning with one loop. Therefore, in the 
identity (|24[) at one loop, the fermion two point func- 
tions on the right hand side have to be the ones at the 
tree level. As a result, the analysis simplifies and gauge 
independence of the fermion mass (pole of the propaga- 
tor) at one loop follows in a simple manner. This is, of 
course, completely in agreement with our explicit eval- 
uation discussed in (j2"0)) . At higher loops, however, the 
analysis is more involved than in covariant gauges. This 
is because in any non covariant gauge involving a distinct 
tensor structure (such as n M ), the structure of S~ 1 (p) be- 
comes quite distinct from that of the tree level two point 
function plj . For example, in two dimensions in the 
axial gauge we can expand the fermion two point func- 
tion at higher loops in a complete basis of Dirac matrices 
T A = (1, 7 m 75 ),/i = 0,1 as 



the covariant gauges where the coefficients have to be 
functions of (m,p M ) which makes the two point function 
to be of the same form as in the tree level theory.) The 
dependence on a new structure makes it much harder to 
analyze the question of gauge independence of the pole. 
As a result, even though we know the Nielsen identity 
to hold at two loop level, we have not yet been able to 
conclude whether the pole of the propagator is gauge 
independent at that order. This is consistent with our 
earlier observation that even though there is no infrared 
divergence at two loops, we are unable to say if the pole 
is gauge independent without explicitly evaluating the 
self-energy at two loops. A general study of the analysis 
of gauge independence directly from the Nielsen identity 
is presently underway. 

To summarize, we have shown explicitly that the pres- 
ence of infrared divergences as well as mass shell singu- 
larities can invalidate the conventional arguments for the 
gauge independence of the pole of the fermion propaga- 
tor. We have shown this explicitly in the case of the one 
loop fermion self-energy in the (two dimensional) mas- 
sive Schwinger model in covariant gauges. We have also 
shown that the leading divergence adds up to zero when 
all the diagrams are summed in the quenched rainbow 
approximation. This suggests that the self-energy may 
be better behaved in a physical gauge. We have shown 
that this is indeed the case in the (space-like) axial gauges 
where the PV (principal value) prescription softens the 
infrared divergence and leads to a gauge independent pole 
of the propagator at one loop. We also find that the in- 
tegrand (in the self-energy) leads to an infrared finite 
integral at two loops although we have not yet been able 
to explicitly evaluate the integral. Therefore, we are not 
yet able to say if the pole of the fermion propagator is 
gauge independent even at two loops. The same conclu- 
sions follow from an analysis of the Nielsen identity up to 
this order and this question is being further studied. If 
the axial gauge does allow a gauge independent definition 
of the dynamical fermion mass in the massive Schwinger 
model, it would indeed be worth trying such an analysis 
in the calculation of the quark mass in QCD. Details of 
our calculations will be published elsewhere. 
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where the coefficients a a are, in general Lorentz covariant 
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